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The secretion of vesicles for intracellular transport often rely on 
the aggregation of specialized membrane-bound proteins into a 
coat able to curve cell membranes. The nucleation and growth of 
a protein coat is a l<lnetlc process that competes with the energy- 
consuming turnover of coat components between the membrane 
and the cytosol. We propose a generic kinetic description of coat 
assembly and the formation of coated vesicles, and discuss Its 
implication to the dynamics of COP vesicles that traffic within the 
Golgi and with the Endoplasmic Reticulum. We show that station- 
ary coats of fixed area emerge from the competition between coat 
growth and the recycling of coat components, in a fashion resem- 
bling the treadmilling of cytoskeletal filaments. We further show 
that the turnover of coat components allows for a highly sensitive 
switching mechanism between a quiescent and a vesicle produc- 
ing membrane, upon a slowing down of the exchange kinetics. 
We claim that the existence of this switching behaviour, also trig- 
gered by factors such as the presence of cargo and variation of 
the membrane mechanical tension, allows for efficient regulation 
of vesicle secretion. We propose a model, supported by differ- 
ent experimental observations, in which vesiculatlon of secretory 
membranes Is impaired by the energy consuming desorption of 
coat proteins, until the presence of cargo or other factors triggers 
a dynamical switch Into a vesicle producing state. 

Transport vesicle | protein Coat | COP vesicles | self- assembly | non- 
equilibrium phase transition 



Introduction 

The plasma membrane and the membrane of cell compartments such 
as the ER and the Golgi continually produce vesicles for cargo trans- 
port. Vesicle formation generally involves specific proteins that ag- 
gregate into semi-rigid coats of dimensions in the lOOnm range, well 
visible by electronic microscopy |T]|2l[3l. The process of vesicle 
formation is now rather well established f4l, and is sketched Fig.l. 
First various cytosolic proteins assemble on the membrane into ele- 
mentary coat-building units, called monomers in the following. The 
membrane-bound monomers then polymerize into coat structures that 
locally bend the membrane and recruit cargo molecules. As the coat 
expands, the coated membrane invaginates until forming a nearly 
spherical vesicle containing cargo |5 , 6 1, that is eventually released 
from the membrane. The coat components soon disassemble and are 
ready to participate to the formation of a new vesicle. 

The coats are classified in three major classes, COPII, COPI and 
Clathrin. Although they involve distinct proteins, the three types of 
coat share many common features, from their size and shape to the 
mechanism by which polymerization, cargo recruitment and mem- 
brane deformation is achieved |4|. Our approach is primarily aimed 
at studying the formation of COPI and COPII vesicles. However, the 
generality and robustness of its outcome suggest relevance for the 
more sophisticated Clathrin coats as well. The assembly of COPs 
and Clathrin, and the fission of COP vesicles can now be reconsti- 
tuted on purified liposomes with a restricted number of components 
EllHllFl- Those experiments point-out the robustness of the coat for- 
mation process. They also confirm that coat polymerization is spon- 



taneous, only driven by weak short range attractions between the 
monomers, while coat disassembly requires the presence of an en- 
ergy source. More precisely, the assembly and disassembly of COP 
coat components follow the cycle of activation - inactivation of a GT- 
Pase protein, Sari for COPII and Arfl for COPI |4|. Once activated, 
the GTPases bind to the membrane and recruit individual coatomer 
complexes (the monomers), that later polymerize into coats 1101 . The 
inactivation of the GTPase, triggered by the hydrolysis of its bound 
GTP, leads to its unbinding from the membrane and to the monomer 
disassembly if the GTPase belongs to a monomer. 

Strikingly, FRAP experiments suggest that the exchange kinet- 
ics of coat components is much faster than the rate of vesicle secre- 
tion nOj . In other words, many futile monomers are released to the 
cytosol during the expansion of a coat. So, while new membrane- 
bound monomers polymerize at the coat periphery, others within the 
coat disassemble and are expelled to the cytosol. Paradoxically, the 
consumption of energy via GTP hydrolysis seems to work against 
coat growth and to prevent vesicles formation. This resembles micro- 
tubules dynamics and by analogy to the treadmilling of microtubules, 
it has been suggested that the competition between growth and un- 
binding may produce stable coats of fixed area I10II1III12III3I . 

In this paper, we investigate theoretically the consequence of fu- 
tile release of coat components on the distribution of size and shape 
of protein coats, and more practically on the amount of secreted vesi- 
cles. In our model (Fig.l), monomers are continuously "dropped" 
onto a membrane and proceed to aggregate into coats of growing 
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Fig. 1 . a) A monomer cycle: activation/membrane binding (Jon), aggregation 
(jp(£)), membrane unbinding via inactivation (joff (^)) or vesiculatlon (jv). b) 
Coat population and the global fluxes in and out of the membrane. 
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size that curve the membrane. Monomers leave the membrane ei- 
ther individually after GTP hydrolysis or collectively as part of a 
completed vesicle. Intuitively, one expects GTPase inactivation to 
decrease the rate of vesicle formation by reducing the lifetime of 
membrane-bound monomers. However, our generic approach reveals 
that a deeper understanding of vesicle secretion requires a quantita- 
tive statistical model. Indeed, we report the existence of a discon- 
tinuous dynamical transition from a quiescent to a vesicle producing 
membrane, upon variation of the rate of GTP hydrolysis. In other 
words, the apparently counter-productive energy consumption that 
favors the unbinding of coat components provides secretory mem- 
branes with a highly sensitive switch to regulate vesicle release, trig- 
gered for instance by a variation of the cargo concentration or the 
mechanical tension of the membrane. 

Description of tlie model 

Our goal is to describe the collective behaviour of a population of 
evolving membrane domains (coats) formed by the aggregation of 
identical units (monomers), which are themselves continuously recy- 
cled between the membrane and a reservoir (the cytosol). Our starting 
point is the course of events depicted in Fig. 1. We consider a patch 
of membrane much larger than the size of individual coats, which is 
subjected to a constant and homogeneous in-flux of monomers Jon. 
The monomers have a finite lifetime on the membrane before being 
recycled to the cytosol, at a rate fcog. While on the membrane, they 
diffuse and eventually aggregate into curved protein coats. Coats that 
manage to reach a critical size leave the membrane as coated vesicles. 

Monomers : The formation of a new monomer on the membrane 
involves a succession of steps (GTPase binding on membrane and ac- 
tivation, and the recruitment of coat proteins) which are not individ- 
ually described in the present model. The rates associated with these 
processes enter a unique parameter, the mean number of monomers 
Jon formed on the membrane per units of time and area. 

Coat growth: Coat expansion proceeds by polymerization of 
monomers at the coat edge. Monomer-monomer binding is spon- 
taneous and results from weak short range interactions. The binding 
energy 7 should be in the range of a few k-^T (k^T is the energy 
available from thermal fluctuations, with fcg the Boltzmann constant 
and T the temperature in Kelvin), since the lOfc^ T provided by GTP 
hydrolysis is sufficient to break the bonds. The polymerization is thus 
thermally reversible and solely driven by the minimization of the free 
energy of the coat. 

Coat structure : Electron microscopy (T]|2l supports the assump- 
tion that the optimal area per monomer sq (~ 100 nm^) and the 
optimal radius of curvature of the coat _Ro (~ 50 nm) are homo- 
geneous within the coat and remain constant during coat growth. We 
thus adopt a model in which the state of a coat is fully characterized 
by a single, slowly varying parameter: the number of polymerized 
monomers £ it contains, taken as a continuous variable for commod- 
ity. All other internal degrees of freedom in the coat (protein density 
and coat shape) are considered to adjust to their optimal configura- 
tion faster than the typical rates of coat growth and GTP-hydrolysis. 
Under these assumptions, a given coat of size £ can be described as 
a spherical cap of constant curvature (defined as the dimensionless 
quantity c = ^so/(47ri?g) ~ 1/20). A full spherical coat (£c^ = 1) 
with these properties contains several hundreds monomers. 

Monomer release : Contrary to the reversible monomer polymer- 
ization, monomer desorption is an energy consuming process driven 
by GTP hydrolysis. It occurs at a rate fcotr, assumed constant here 
for simplicity (see Supporting Information (SI) for a discussion of 
this assumption). This rate can be estimated from FRAP experiments 
(T0l[T2l[Tl[T5l, feoff ~ 0.1 - 10 s~\ Under the assumption of con- 



stant protein density in the coat, the dissociation of monomer from 
the coat is followed by a rapid rearrangement of the coat structure 
and by a slight shrinkage of the coat; monomer inactivation thus op- 
poses coat expansion. 

Vesicle release : A mature coat containing a number £v {£v(? ^ 
1) of bound monomers forms a nearly closed sphere connected to the 
rest of the membrane by a thin neck. At this stage, the coat cannot 
grow further, and is eventually released into the cytosol as a coated 
vesicle. The details of the scission mechanism vary between classes 
of coat, and may involve additional proteins (|4|). Here, we merely 
assume that once a domain reachs the critical size £ = £v,\\. leaves 
the membrane as vesicles at a constant rate fcv 

Membrane properties: The curvature of the coat imposes a de- 
formation to the membrane which is opposed by membrane tension 
(16!|. The membrane tension a thus favors coat depolymerization, 
and can have a sizable effect on coat growth if it is larger than 
j/isoV^v) — 10~^ J/m^ (A detailed model for the elastic proper- 
ties of the protein-covered membrane is discussed in SI, section II.B). 
Tensions of the Golgi, the ER, and the plasma membranes are typical 
in the range a ~ 10^'' — 10^* J/m^ 1 171 . and may thus play a role 
in vesicle secretion. Hereafter, membrane tension will be expressed 
in natural units: a = aso ~ 10~^ — Ik^T. 

The population of membrane coats is characterized by its size 
distribution n{£) . The mean concentration of coats of size £ (between 
£ and £ + d£) at a time t is n{£, t)d£, and the mean concentration of 
isolated monomers is ni(f). Our purpose is thus to compute, n{£) 
and rii at steady state, for given values of the parameters 7, a, feofr, 
Jon, fcv and £v More practically, we will compare the fluxes of coat 
elements leaving the membrane as inactive monomers Joff and as 
part of a vesicle Jv (Fig. 1). 

Thieoretical frameworic 

Monomer fluxes and conservation relations. In this section, 
we derive the kinetic equations for the evolution of the coat size dis- 
tribution n{£). The monomer cycle can be divided into four steps, to 
which correspond four different fluxes, as shown Fig. I. 

- Jon (i) is the in-flux of single monomers binding to the membrane, 
taken as an input in our model. 

- jp{£, t) is the flux of monomer joining domains of size £, a bal- 
ance between polymerization and depolymerization for this domain 
size. Integrated over the entire population, it gives the total flux of 
monomers incorporated into domains Jp = /, ^ dijp{£,t). 

- joff {£, i) is the flux of monomers expelled from domains of size 
£ into the cytosol after GTP hydrolysis. Under the assumption of 
uniform release, it is given by joff (^) = koffn(£)£. Integrated over 
the entire population, it gives the total flux of individual inactive 
monomers leaving the membrane Joff = L^ d£ joff {£,t). 

- jv (t) is the flux of mature coats released as vesicles. Introducing the 
rate of vesicle formation fev, we have jv = fcvn(£v). The total flux of 
monomers leaving the membrane as part of a vesicle is Jv (i) — jv£v 

The evolution of the coat size distribution satisfies (see SI): 

dtn{£) = ^di(M£)^joff{£)), [1] 

and that the polymerization current jp reads: 

JpW = -kpni{den{£) +n{£)deAE{l)). [2] 

jp is proportional to the density of available monomers ni, and to 
the rate of monomer binding onto coats fep (see below). It contains 
a diffusive term accounting for random polymerizations and depoly- 
merization induced by thermal noise, and a convective term describ- 
ing the drift of monomers toward domains of lower energy, driven by 
the "force" —diAE. The free energy difference AE (in fc^T units) 
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between a coat of size I and I isolated monomers diffusing on the 
membrane is obtained treating the coat as a rigid spherical cap 1 18| 
(see also SI): 



where /i = In ni + 7\/l — c'^ + <^(? 



[3] 
[4] 



is a chemical potential including the entropy of the freely diffusing 
monomers, see SI. 

The net current j{£) — jp{£) — jos{£) accounts for polymeriza- 
tion and desorption. It can be written in terms of an effective energy 

Ei£) = AE{1) + ^^{l^ - 1) : 



j = —kpni{den + ndiE), 



with 



E{£) = j^£{l-c'^e) + T,{ni)e'^ - ^(m) £ + const 

S(ni 



[5] 



[6] 



- 2 fcoff 1 ,-_, 

(^c +— . [7] 

ifcp ni 



This equation introduces an effective tension E that illustrates the 
fact that desorption of inactive monomers and membrane tension for- 
mally play the same role in hindering coat maturation and vesicle 
secretion. Note that more generally, the coatomer binding or inac- 
tivation rates may depend on the coat size, in which case monomer 
desorption enters the effective energy as J d££{koB{£)/kp{£)). 

Finally, the monomer influx Jon and the flux of secreted vesicle 
jv are accounted for via the boundary conditions (see SI) 



Jp{^) — '-'on Jp, 



[8] 
[9] 



Steady state. At steady state, all fluxes are balanced and dtn — 0. 
Eas.( IU9t reduce to two conditions to be satisfied by ni and n{£): 



j{^) ~ h = constant. 

t^on — t-'off ~r Ov , 



[10] 

[11] 



The former equation enforces that the size distribution is constant, 
and the latter that the flux of monomer binding to the membrane bal- 
ances the flux of monomer leaving the membrane, either after inacti- 
vation or by vesiculation. 

Results 

In this section, we focus on the steady state of a membrane receiv- 
ing a constant in-flux of monomer, each having a finite lifetime at 
the membrane. The full characterization of the coat population and 
of vesicle secretion follows two steps. First, the stationary distribu- 
tion of coat size n{£) is computed for a given concentration of free 
monomers n\ with Eq.lllOt. Second, n-i is self-consistently derived 
for a given monomer in-flux Jon by imposing that the in-flux matches 
the total monomer out-flux (Eg. d lit). While the first step relies en- 
tirely on the properties of the free energy landscape E{£), the second 
introduces collective effects emerging from the competitive growth of 
many domains, which ultimately give rise to the "secretory switch". 

Effective energy landscape and steady-state distribution. 

Apart from thermal fluctuations, a coat is driven toward growth or 
shrinkage by the effective "force" —diE[£) (Eq.l[5]l). Coat growth 
is thus formally analogue to thermal diffusion along the effective en- 
ergy landscape E{£) (Eq.l|6]l). The analytical expression of the coat 



size distribution is given in the SI. Several different regimes can be 
distinguished under increasing monomer concentration (Fig.2). 

The energy landscape illustrates the interplay between antago- 
nistic effects; short-range attractions between monomers promotes 
polymerization, while the entropy of the free monomer favors their 
dispersion. The 7 and ^ terms in Eq.([6f reflect this competition. 
Furthermore, monomer inactivation and unbinding, and membrane 
mechanical tension hinder coat growth (kos and a respectively, com- 
bined in the effective tension E, Eq.O). Depending on the monomer 
concentration, the landscape may show a local maximum at a small 
size £-a (Fig.2b-d), which indicates a nucleation process. Coats must 
reach the critical size £-a (through fluctuations in the pool of free 
monomers) in order to consistently grow further, and the rate of nu- 
cleation is controlled by the height of the energy barrier. The effec- 
tive energy may also show an barrier to vesiculation for large coat 
size £ ~ £v (Fig.2b-c), indicating that large coats are suppressed by 
the effective membrane tension E. A local minimum then exist for 
an intermediate size f , corresponding to kinetically stable coats. 
Low monomer concentration. At low monomer concentration, the 
large energy barrier to vesiculation at £ = ^v prevents coats to mature 
into vesicles ( Jv — 0) (Fig.2, states a and b). The coat size distribu- 
tion resembles a distribution at thermal equilibrium: n{£) ~ e^^'*'. 
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Fig. 2. Effective energy landscape E(£) (left column, in k^T units, from 
Eq.lU) and the corresponding coat size distribution n{(.) (right column, in n'^^ 
units, from Egs. lSllOl ), for different values of the free monomer density n\. 
ni/nf^" = 0.986 (a - first row green curve), 1.039 (b - first row red), 1.054 
(c- second row) and 1.160 (d-third row). The variation of E{£) and n{£) upon 
a slight increase of ni above the given value are shown in blue dashed lines. 
Other parameters are 7 = 5 k-^T, a = 0, k^ = O.lSs^-'^, fcoff = 0.125^-'^ 
and, £v = 500. With those values, nf^" = 0.009. 
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and the membrane follows a classical scheme common to many self- 
assembling systems (e.g surfactants in solution U?]). The local en- 
ergy minimum at £* appears above a critical concentration rii™'^, ana- 
logue to the "critical micellar concentration", or "cmc" at which sur- 
factants in solution start forming micellar aggregates (see 1 19] and 
the SI). For ni < n™'^ (Fig. 2a), entropy dominates and the effective 
energy increases monotonously with the coat size £. Monomer ag- 
gregation is unfavorable, and the membrane contains mainly single 
monomers and few small transient domains formed by fluctuation. 
For ni > n""'^ (Fig. 2b), long-lived coats can nucleate, at a rate 
fixed by the nucleation barrier, and grow up to the optimal size I*. 
Maturation into coated vesicles (i = iv) is prevented by monomer 
desorption and membrane tension. 

Larger monomer concentration. The height of the energy barrier to 
vesiculation at £ — l^ decreases with increasing monomer concen- 
tration. When it falls below the nucleation energy barrier (Fig.2c), 
domains may mature into fully-formed vesicles and vesicle secretion 
becomes increasingly probable. The optimal coat size £* increases 
with ni, and eventually exceeds the critical size for vesiculation £v 
at high monomer concentration (Fig. 2d), under which conditions any 
nucleated domain matures into a fully formed vesicle. 

The growth of individual coats is controlled by the amount of 
free, active monomers ni. On the other hand, the pool of free 
monomer is depleted by their binding onto growing coats and is thus 
influenced by the coat population. As we shall see next, this feed- 
back induces remarkable collective effects within the coat population, 
which presents a discontinuous transition between a state of arrested 
growth and a state of abundant vesiculation within a narrow range of 
kinetic parameters. 



Vesicle secretion is controlled by collective effects. The so- 
lution of the coupled Eas.l ilOliTt is graphically represented on Fig. 3 
as the intersection of the monomer in-flux Jon and total out-flux 
JoS + Jv. It may fall in four different regimes (a to d), correspond- 
ing to the four distributions plotted in Fig. 2. In a wide range of pa- 
rameters (see below), Job displays the remarkable property of being 
non-monotonous, with a sharp peak at a critical concentration of free 
monomers. This behavior dramatically influences the membrane's 
ability to secrete vesicles. Indeed, a given monomer in-flux may cor- 
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Fig. 3. The fluxes of monomers leaving the membrane at steady state, as a 
function of the density of active monomer ni (in n™"^ unit). Jv (blue) is the vesic- 
ulation flux and Joft (red) the sum of vesiculation and inactivation. At stationary 
state, the total outgoing flux balances the incoming flux ( Joff + Jv = Jon) (hori- 
zontal dashed line). Fluxes are in fcp(n5;™'^)^ unit and the blacl^dots correspond 
to the three different states (b, c, d) depicted in Fig.2. 



respond to three distinct dynamical states of the membrane. We will 
show below that regimes b and d represent respectively a quiescent 
membrane and a membrane secreting large amount of vesicles, while 
regime c is dynamically unstable. The secretory membrane thus con- 
stitutes a bistable dynamical system able to abruptly switch vesicle 
secretion on and off at prescribed monomer turnover rates. 

Since all membrane-bound monomers are inactivated with the 
same rate kas, the total flux of monomer leaving the membrane af- 
ter inactivation Joff is directly proportional to the total amount of 
monomer on the membrane. The peak of Joff in Fig. 3 stems from the 
complex relationship between the concentration of isolated monomer 
n\ and the total amount of coat components on the membrane. 
No vesicle secretion regimes a and b. If no coat can form (low 
monomer concentration: state a), the monomers out-flux is domi- 
nated by the desorption of free monomer: Joff — fcoffni. If coats can 
form, but do not mature into vesicles (state b), the out-flux is domi- 
nated by the desorption of monomers belonging to coats of size £* : 
Joff ~ koft£*n{£*). In this regime, monomers reaching the mem- 
brane tend to join a coat and the density of free monomer is almost 
insensitive to the fluxes: (ni ~ n{£*y^' with £* 2> 1, see SI). 
A small increase of ni requires a pronounced increase of the to- 
tal amount of coat material on the membrane, which explains the 
sharp rise of the total monomer out-flux Joff with ni in Fig. 3. In 
this regime, the optimal coat size is also insensitive to the monomer 
fluxes, and is obtained from the minimization of the effective energy 
E (Eq.H) 

(* « ( 7T7-^ ) , [12] 

VE(n5;""<=)y ' ^ ^ 

Tlie unstable regime c. For intermediate monomer density, 
metastable coated pits have a high probability to grow into fully 
formed vesicle owing to the small barrier to vesiculation (Fig. 2c). 
The rate of vesicle formation increases with n\, so the total amount 
of membrane-bound material actually decreases with increasing con- 
centration of free monomer. This is shown by the dashed blue line in 
Fig.2c, and explain the decrease of Joff in Fig. 3. This situation can- 
not be maintained at steady state, and spontaneously evolves toward 
either state b or d. 

Steady vesicle secretion regime d. If the monomer concentration on 
the membrane is large, the effective coat energy exhibits a nucleation 
barrier but no intermediate minimum (Fig. 2d). After nucleation, a 
coat grows at nearly constant velocity until reaching the critical size 
£v where it remains trapped for a time 1/fcv before being released as 
a vesicle. In this regime, both Joff and Jv increase with ni, Fig.3. 

The bistability exhibited by the coats dynamics relies on the ex- 
istence of an unstable steady state and holds as long as there ex- 
ist a (meta)stable coat of intermediate size. This feature is con- 
served even if the effective energy contains higher order terms, to 
be expected if the ratio of monomer dissociation to binding rates 
(fcoff /ftpTii) increases with the coat size (see SI). Furthermore, the 
switch exists if there is a metastable-state within the accessible size- 
range: £* /£v{~ £*c?) < 1. From Ea.J12t. this condition amounts 
to E > 7c^(~ 10~^). The effective tension E (Eq.Q) accounts 



both for the membrane mechanical tension (ctc 



10" 



10^-') 



and the ratio (fcoff/fcpni). The binding rate is assumed to be lim- 
ited by monomer diffusion, and is expected to be of order the inverse 
monomer diffusion time over its own size (fcp ~ D/ sq ~ 10^s~^, 
with D ~ pm} I s the membrane diffusion coefficient). With a dis- 
sociation rate fcoff — ls~^, and a monomer density n\ ~ 1% (or 
fcpWi ~ 10^s~^), we find that secretory membranes are well into the 
bistable regime (E ~ 10~^ ^ 1(?), and should exhibit the secretory 
switch discussed below. 
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Discussion 

The growth of coated pits and the secretion of coated vesicles result 
from a kinetic balance between the polymerization and the inactiva- 
tion of coat components. It is thus to be expected that coat maturation 
can only proceed if the coatomers turnover at the membrane is suffi- 
ciently slow llOllll|[T3l . Our model goes beyond this intuitive anal- 
ysis, and shows that secretory membranes are able to abruptly switch 
between a quiescent and a vesicle producing state upon a slowing 
down of coatomer recycling. Switch-like behaviors are clearly ad- 
vantageous for biological systems. The highly non-linear nature of 
a switch confers evident robustness with respect to the noisy envi- 
ronment, and allows for a precise regulation of the system's activity. 
The biological consequences of the "secretory switch" are discussed 
below, together with the influence of important factors, such as the 
density of cargo or the membrane tension, in regulating the activity of 
secretory membranes. We also show how some apparently unrelated 
observations on COPs vesicles naturally fit into our global picture of 
coated vesicle secretion. 

The "secretory switch". Consider a secretory membrane receiv- 
ing a fixed amount of coatomer per unit time. As coatomers accumu- 
late and aggregate, the membrane eventually reaches a steady state 
in which the flux of coatomer leaving the membrane (by inactivation 
or vesicle secretion) balances the in-flux. If the in-flux is low, the 
membrane is covered by monomers and stationary coated pits, the 
latter being prevented to mature into fully-formed coated vesicles by 
the combined effect of coatomer recycling and the energy associated 
to membrane deformation. As the in-flux increases, the coated pits 
to monomer ratio increases (Fig.2a, dashed blue line), but as long as 
the membrane remains in the "stationary pits" regime (b - Fig. 3) the 
size and shape of the coated pits is weakly sensitive to the in-flux and 
no vesicle is produced. However, beyond a threshold value of the in- 
flux, state b disappears (Fig.3). Coated pits are not kinetically stable, 
and after a fast transient regime that sees the release of the previously 
stable coated pits, the membrane settles into a state of steady vesicle 
secretion (state d). 

If the in-flux is now reduced, the membrane remains in the se- 
cretory state d, which possesses a stable branch for smaller in-flux 
(Fig.3). The secretory regime disappears below yet another critical 
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Fig. 4. number of secreted vesicle per ^im? .s (= jv x 1/^m/so) as a func- 
tion of the rate of coatomer desorption from the membrane feoff. The transition 
to vesiculation is discontinuous, and is characterized by an hysteretic cycle (ar- 
rows), with high and low turnover thresholds. The parameters are the same as 
those used for Fig.2 and 3. 



threshold and the systems jumps back into the quiescent state b. 

The coat population can thus undergo discontinuous dynamical 
transitions, characterized by an hysteretic cycle, between the two 
non-equilibrium steady states b and d. In other words, the secretory 
membrane works in an all-or-nothing fashion and can switch vesicle 
production on and off within a narrow range of control parameters 
such as the GTPase activation and inactivation rates, and membrane 
tension. Fig.4 shows the rate of vesicle secretion as a function of 
the GTP hydrolysis rate fcoff . The discontinuous transition between 
quiescent and vesicle-producing membrane is clearly apparent, and 
is characterized by two hydrolysis rates (high and low thresholds). 
Between these two rates, the secretory membrane may be in either 
state, depending on the system's history (hysteresis). 

Regulation of vesicle secretion by cargo. The adsorption flux 
Jon, and the desorption rate fcoff, of coat components at the mem- 
brane tightly control vesicle secretion. Recent fluorescence experi- 
ments on COPs coat suggest that these rates vary with the amount 
of cargo present at the membrane. For COPI, the presence of extra 
cargo leads to significant increase of the amount of coat components 
at the Golgi membrane, which could reflect either the increase of Jon 
or the decrease of fcofr 1121 . For COPII, FRAP experiments show 
that the coatomer exchange rate between the ER membrane and the 
cytosol is doubled in the absence of cargo. This has been attributed 
to the increase of fcoff with decreasing cargo density |T5 1. 

Our model predicts that vesicles can only be secreted if the re- 
cycling rate kos is below a critical value (Fig.4). By increasing the 
lifetime of the coatomers at the membrane, the presence of cargo is 
thus expected to promote vesicle secretion, and a minimal amount of 
cargo at the membrane might actually be required for transport vesi- 
cle to be secreted. The two, high and low, recycling thresholds of 
Fig.4 would then corresponds to two critical cargo densities (low and 
high, respectively). 

Considering that newly synthesized cargo is brought to the 
membrane at a (slow) steady rate and is removed by vesicula- 
tion, membrane-bound cargo accumulates in the no-secretion regime, 
thereby decreasing fcoff and moving the system toward the secretion 
regime. Above the high cargo density threshold, the coat-machinery 
abruptly escapes the stationary pits regime and switches to vesicle 
production (Fig.4). This results in a decrease of membrane-bound 
cargo, which increases koS and moves the system toward the quies- 
cent state. Below the low cargo density threshold, vesicle secretion 
is switched off, letting the cargo accumulate until the high-density 
threshold is reached and vesicle production is resumed, starting a 
new cycle. Under constant cargo in-flux, the system should thus pe- 
riodically switch between quiescent phases and phases of vesicle se- 
cretion, following the hysteretic loop of Fig.4. If cargo synthesis is 
irregular, the membrane waits for sufficient accumulation of cargo 
between transient residences in the secreting state, where the accu- 
mulated cargo is released. The vesicle flux over time should then ap- 
pear as an irregular pulsed signal. Recording the total vesicle out-flux 
over an extended patch of secretory membrane over time should thus 
be of high interest. Oscillatory or pulsed vesicle secretion would be a 
signature of the secretory switcli uncovered by our theoretical analy- 
sis (within our framework, steady secretion would indicate that cargo 
synthesis is sufficiently fast to compensate the secreted cargo). 

This accumulator mechanism would provide functional effi- 
ciency to the secretory membrane, as it would prevent the futile de- 
livery of empty vesicles. Strikingly, analysis of COPI vesicles in 
mutant cells where arf I is unable to hydrolyse GTP have revealed a 
much lower cargo content than in normal cells 1201 . This observation 
supports our prediction. Indeed, in the absence of GTP hydrolysis 
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(fcoff = 0), the switch is gone and vesicle secretion remains "on" 
regardless of the available amount of cargo. 

The existence of an hysteretic cycle can also positively impact 
on the rate of cargo delivery. At steady state, the flux of cargo deliv- 
ered by a secretory membrane that would not possesses an unstable 
regime would automatically adjust to the flux of synthesized cargo. 
Here, and while the system is traveling along the secretion branch of 
the hysteretic cycle, the flux of secreted cargo is mainly controlled 
by the kinetics of coat formation, and can potentially be much larger 
than the rate of cargo synthesis. 

Effects of the membrane tension: Regulation of vesicle 
formation and coat flattening. Our calculation has shown that 
membrane tension plays essentially the same role as coatomer re- 
cycling in opposing vesicle secretion (Eq.l[6l). Vesiculation is only 
possible below a high tension threshold, and the oscillatory behavior 
described above may also result from variations of tension. Secre- 
tion removes membrane area from the organelle and may increase its 
tension, while the fusion of incoming vesicles 1211 or other regula- 
tory mechanisms 1181 , dynamically relaxe tension. Vesicle secretion 
would thus occur only when enough membrane area has been accu- 
mulated to relax the tension. Such a mechanism suggests a coordi- 
nation of purely mechanical origin between absorption and release 
of vesicles, and could prevent the uncontroled shrinking of secretory 
compartments. Remarkably, the Golgi strikingly crumbles in cells 
where the GTP hydrolysis in COPI-coat is rendered inoperant |20|. 

In Eq.llS), the effect of membrane tension was computed assum- 
ing that the coat rigidity k (unit of energy) was sufficiently strong 
to impose the domain curvature, which remained constant regardless 
of the coat size and the membrane properties. This simplification 
ceases to be valid under high membrane tension (a ~ h/Rq, where 
Ro is the radius of curvature of a tensionless coat). Higher tensions 
result in a flattening of the coat (see SI, section II.B). Beyond a ten- 
sion threshold a > k/{2RI) (~ 10~*J/m^ for the coat rigidity 
K = lOOfcg T), the formation of a closed sphere becomes impossible 
and protein aggregates should grow as flat patches. High membrane 
tension would thus favor the formation of flat coatomer aggregates 
which size is limited by coatomer recycling |22|. Such "flat lattice" 



are indeed observed for Clathrin coats at the basal membrane of ad- 
hered cells [231, where adhesive proteins are expected to generate 
high membrane tensions. 

Concluding remarks. The model presented here is the simplest 
implementation of a kinetic model of coated vesicle formation where 
coat growth competes with the inactivation of coat components. The 
secretory switch revealed by our work is very robust and relies solely 
on the (non-equilibrium) Thermodynamics of coat formation. Be- 
yond the qualitative agreement with experimental findings discussed 
in the previous section, we hope that future experiments can fur- 
ther test some of our predictions, which include: ;) the existence of 
metastable domains of intermediate size, //) the role of membrane 
tension in preventing the formation of curved protein coat, and even- 
tually its involvement in the formation of flat lattices. Hi) the oscilla- 
tory or pulsed secretion of vesicles in time. 

In this study, we have used the crudest possible description of the 
coat structure, thus avoiding to deal with structural details of specific 
protein coats. Further experimental observation, e.g on biomimetic 
systems, could motivate the building of models focusing on the dy- 
namics of a single coat. Supplementary degrees of freedom for the 
coat shape could be considered, allowing for the competitive growth 
of structures of various morphologies (tubules L6J, spherical caps and 
flat lattices 1231 ). observed in living cells and biomimetic systems. 

In the same spirit, the heterogeneties of the coat structures may 
be included in the model. The coupling beween the GTP hydrol- 
ysis rate and the curvature (COPI) or the degree of polymeriza- 
tion (COPII) revealed in recent experiments I14II24I . suggest that 
monomers inactivation may be enhanced at the coat center and lim- 
ited at the boundary. One could then imagine the formation of layer 
of active monomers preventing coat disassembly 1141 1241 . Such 
properties suggest a strong analogy with microtubules | Hi, and it is 
tempting to imagine exotic growth dynamics with shrinking cascades 
such as those observed for microtubules 1251 . 
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